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We use the electronic spin of a single Nitrogen- Vacancy (NV) defect in diamond to observe the 
real-time evolution of neighboring single nuclear spins under ambient conditions. Using a diamond 
sample with a natural abundance of ^^C isotopes, we first demonstrate high fidelity initialization 
and single-shot readout of an individual ^'^C nuclear spin. By including the intrinsic nuclear 
spin of the NV defect in the quantum register, we then report the simultaneous observation of 
quantum jumps linked to both nuclear spin species, providing an efficient initialization of the two 
qubits. These results open up new avenues for diamond-based quantum information processing 
including active feedback in quantum error correction protocols and tests of quantum correlations 
with solid-state single spins at room temperature. 
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Nuclear spins are attractive candidates for solid-state 
quantum information storage and processing owing to 
their extremely long coherence time^^^. However, since 
this appealing property results from a high level of isola- 
tion from the environment, it remains a challenging task 
to polarize, manipulate and readout with high fidelity 
individual nuclear spinJ^. A promising approach to over- 
come this limitation consists in utilizing an ancillary sin- 
gle electronic spin to detect and control remote nuclear 
spins coupled by hyperfine interactiorPtiSI, this con- 
text, the NV defect in diamond has recently attracted 
considerable interest because its electronic spin can be 
polarized, coherently manipulated, and readout by opti- 
cal means with long coherence times, even under ambient 
conditionJli'. The NV's electronic spin thus behaves as 
an ultrasensitive magnetometer at the nanoscal^', pro- 
viding a robust interface to detect and control nearby 
nuclear spins in the diamond lattice. This approach has 
been used in the past years to study the coherent dynam- 
ics of mu lti-sp in systemJ^^, to perform universal quan- 
tum gates^I^nSI and to develop few-qubits quantum reg- 
isters, where shigle nuclear spins are used as quantum 
memorie d^l^^l^'^ i A second-long coherence time was re- 
cently demonstrated for a single ^'^C nuclear spin weakly 
coupled to a single NV defect in an isotopically puri- 
fied diamond sampl^i2]_ This result, combined with the 
ability to perform spin-photon entanglement^^ and two- 
photon interf erence from distant NV defects at low tem- 
peratur j^^ l ^^ l, makes single spins in diamond a promising 
building block for quantum repeaters and long-distance 
quantum communications. 

However, advanced quantum algorithms such as quan- 
tum error correction protocols require high fidelity ini- 
tialization and single-shot readout over multiple qubita^ii. 
Alo ng the line of recent works directed towards this 
gpa,]SHini22]^ we first report high fidelity single-shot read- 
out of an individual ^^C nuclear spin by using the elec- 
tronic spin of a single NV defect as an ancillary qubit in 



a diamond sample with a natural abundance of ^ C iso- 
topes (1.1%). Repetitive readout indicates a polarization 
lifetime exceeding seconds at moderate magnetic fields, 
which illustrates the robustness of the -'^^C nuclear spin 
state. Then, wc demonstrate efficient initialization of two 
nuclear spin qubits in a well-defined state by adding the 
intrinsic ^^N nuclear spin of the NV defect in the quan- 
tum register. 

The spin system considered in this study is depicted 
in Figs, [ija) and (b). The electronic spin {S = 1) of 
a single NV defect is coupled by hyperfine interaction 
with both its intrinsic "'^^N nuclear spin (/ = 1) and a 
neighboring ^^C nuclear spin (/ = 1/2). A permanent 
magnet placed on a three-axis translation stage is used 
to apply a static magnetic field with controlled amplitude 
along the NV defect axis and the spin transition between 
the fUs = and nis = — 1 electron spin manifolds is co- 
herently driven through microwave (MW) excitation. As 
shown in Fig.jljc), the hyperfine structure of the s pin sys- 
tem, recorded through pulsed-ESR spectroscopy's^ ex- 
hibits six nuclear-spin conserving transitions (see Sup- 
plementary Information). From this spectrum, recorded 
for a magnetic field magnitude B — 2000 G, we extract 
the projected strength of the ^^C hyperfine interaction 
= -^zz = 258 ± 10 kHz. We obtain further qualitative 
information of the hyperfine interaction through dynamic 
polarization measurements at the excited-state level anti- 
crossing (ESLAC), while applying a static magnetic field 
near 510 G along the NV axi^SlT The ^'^C polarization 
efficiency was estimated by using the Fourier transform 
of the free-induction decay (FID) signal measured by ap- 
plying a Ramsey sequence ^ — t — ^ to the NV defect 
electron spin. Figure[ljd)) shows the FID signal recorded 
at the ESLAC. Since the ^^N nuclear spin is perfectly po- 
larized, the characteristic beating is linked to the weakly 
coupled ^^C nuclear spin. The Fourier transform of the 
FID signal indicates a relatively high polarization effi- 
ciency = 40 ± 10%, which suggests that the ^'^C quan- 
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FIG. 1: Spin system, (a), Energy-level diagram of the NV 
defect as a function of the strength of a static magnetic field 
B applied along the NV defect axis. Level anti-crossings in 
the ground state (GSLAC) and in the excited state (ESLAC) 
are highlighted, (b), Hyperfine structure of the rus = and 
rUa = —1 electron spin manifolds for a NV defect coupled 
with its intrinsic ^^N nuclear spin (nuclear spin projection 
mj(N)) and with a nearby single ^"'C nuclear spin (nuclear 
spin projection f or ^). The hyperfine sublevels are denoted 
as |mj(N),ti) Bind the blue arrows indicate the microwave 
(MW) transitions used for single-shot readout measurements, 
(c), Optically detected ESR spectrum recorded for a mag- 
netic field B ~ 2000 G. The ^^N hyperfine interaction leads 
to a splitting of 2.16 IVIHz between ESR frequencies associ- 
ated with different ^''N nuclear spin projections. These lines 
are further split by A^z = 258 ± 10 kHz through hyperfine 
coupling with a nearby ^^C. (d), FID signal of the NV defect 
electron spin recorded at the ESLAC {B ~ 510 G) showing a 
coherence time T2 = 2.9 ± 0.1 fis, in the range expected for a 
diamond sample with a natural abundance of '^^C. (e)-Fourier 
transform of the FID signal showing a significant polarization 
(~ 40%) of the ^^C nuclear spin in state jf). 



tization axis is close to the NV defect axis in both the 
ground and excited state^^U, The anisotropic component 
of the hyperfine tensor Aani is therefore assumed to be 



much smaller than Azz- Since the ^^C nuclear spin gets 
polarized in jf), this measurement also provides unam- 
biguous identification of ea ch ES R frequency to a given 
nuclear spin state, |t) or IDI^MII. 

In the spirit of previous work^^Es!^ projective single- 
shot detection of the ^'^C nuclear spin state is achieved 
by accumulating the NV defect photoluminescence (PL) 
while repeating the sequence depicted in Fig. [2]Ja). The 
NV defect electron spin is first initialized into the — 
sublevel through optical pumping. A controlled not 
(CNOT) gate is then applied to induce an electron spin- 
flip conditioned on the ^"^C nuclear spin state. Finally, 
the resulting electronic spin state is optically readout by 
applying a 300-ns laser pulse. This sequence is repeated 
many times in order to increase the signal to noise ratio. 
The CNOT gate is experimentally realized by applying 
narrowband IvIW 7r-pulses on the electronic spin, which 
selectively drive the ESR transition for a given -'^^C nu- 
clear spin state, e.g. \],). In order to take advantage 
of the full ESR contrast and to get rid off any quan- 
tum jumps linked to the ^''N nuclear spin^, three MW 
sources are used for driving simultaneously "'^'^C nuclear 
spin state-selective transitions from each hyperfine sub- 
levels linked to the ^''N nucleus (Figs.[l]^b) and (c)). 

A typical PL time trace recorded while continuously re- 
peating the sequence is shown in Fig. ^h) . For each data 
point, the PL signal is accumulated during Tf, = 120 ms, 
corresponding to approximatively 20000 repetitions of 
the readout sequence. The signal exhibits well-defined 
quantum jumps linked to the evolution of the ^'^C nu- 
clear spin state. Indeed, when the nuclear spin is in state 
the CNOT gate flips the NV defect electron spin, 
TTig — > rus = —1, and a low PL signal is observed 
{dark state) owing to spin-dependent PL of the NV de- 
fect. Conversely, when the nuclear spin is in state |t) , the 
electron spin remains in the = sublevel at each rep- 
etition of the sequence and a high PL signal is observed 
{bright state). Nuclear spin flips are therefore evidenced 
in real-time as abrupt jumps between two distinct val- 
ues of the PL signal. For a magnetic field B = 1610 G 
applied along the NV defect axis, we infer the charac- 
teristic relaxation times of the ^^C nuclear spin while 
applying repetitive readout Ti^(brig/it) = 2.4 ± 0.1 s and 
Tl,l(dark) = 1.5 ±0.2 s. 

To estimate the readout fidelity, the ^^C nuclear spin is 
first deterministically initialized in a given state through 
a single-shot readout measurement. By introducing an 
initialization threshold A/i,^ (resp. A/i,^), photon count- 
ing events such that A/" < A/i,^ (resp. A/" > A/i_^) 
are post-selected, corresponding to an initialization in 
state 11) (resp. |t)). Using A/i,^ = 615 counts and 
Afi,^ = 845 counts, the initialization fidelity exceeds 99% 
for both nuclear spin states (see Fig.[2]jd)). We note that 
in principle the initialization fidelity can be chosen arbi- 
trarily high by decreasing (resp. increasing) A/i,^ (resp. 
A/i,^), at the price of a high number of lost events. After 
successful initialization, a subsequent readout measure- 
ment is performed allowing to build the nuclear-spin de- 
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where Afr is the readout threshold. For Nr,opt = 
735 counts, corresponding to the maximum overlap be- 
tween the two photon-counting distributions, we extract 
J"^ = J"^ = 96 ± 1.2% (Fig. Ijd)). This fidelity could 
be significantly improved by increasing the collection ef- 
ficiency with diamond photonic nanostructures^® . In ad- 
dition, the selective MW 7r-pulses used for the CNOT 
gate have a duration of 4 /is, corresponding to a spectral 
width of 130 kHz. Given the inhomogeneous linewidth 
of the ESR signal combined with the hyperfine coupling 
strength, the 7r-pulses are therefore not perfectly selective 
leading to a decreased contrast of the projective measure- 
ment, which degrades the readout fidelity. This limita- 
tion could be overcome by using a CVD-grown diamond 
sample isotopically enriched with ^^C atoms, in which the 
inhomogeneous dephas ing ra te of the NV defect electron 
spin can reach few 

We now study the processes leading to nuclear spin 
depolarization. For a magnetic field B applied along the 
NV defect axis (z), the ground-state spin Hamiltonian 
reads T-L = + InBIz + S • • I, where includes 
both the pure electronic spin terms and the hyperfine 
interaction with the intrinsic ^"^N nuclear spin, 7„ is the 
^^C gyromagnetic ratio and A its hyperfine tensor. In the 
secular approximation, this Hamiltonian simplifies to 



FIG. 2: Single-shot readout of a single ^^C nuclear spin, (a). 
Logic diagram of the single-shot readout scheme and corre- 
sponding experimental sequence. For all experiments, the du- 
ration of the TT-pulses is set to 4 /xs. The 300-ns laser pulse is 
used both for spin-state read-out and to achieve an efficient 
preparation of the NV defect electron spin in the rris — 
sub level (|0)e) at each repetition of the sequence, (b), PL 
time trace showing quantum jumps of the ^^C nuclear spin 
state. The solid line is a fit with a two states hidden Markov 
model from which the relaxation time Ti of the nuclear spin 
state is extracted, (c), Normalized nuclear-spin dependent 
photon counting distributions 5^(4.). The solid lines are data 
fitting with the formula given in the Supplementary Informa- 
tion, (d). Single-shot readout fidelity as a function of 
the readout threshold. The initialization thresholds A/i4,A/i,-|- 
and the optimized discrimination threshold Nr,opt are indi- 
cated with dashed lines. The solid lines are extracted from 
the fits in (c). We note that the initialization fidelity in state 
l;) (resp. It)) is given by Ti{M^,i) (resp. 



pendent photon counting distributions As shown 

in Fig. [2]^c), the distributions linked to each nuclear spin 
state can be clearly distinguished and the readout fideli- 



ties -^■f(^) are defined as 



(1) 
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where Aani = (^L + -^lyV^^^ tan(/i = Azy/Azx and 
I± = Ix ^ ily The anisotropic component Aani of the 
hyperfine tensor therefore induces nuclear spin flips, lead- 
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FIG. 3: Nuclear spin relaxation time T\ as a function of the 
strength of a magnetic field applied along the NV defect axis 
for the bright state (red) and the dark state (blue) . The solid 
lines are data fitting with a simple model including nuclear 
spin flips induced by the transverse component of the hyper- 
fine tensor and electron-nuclear spin fiip-fiops at the GSLAC 
and ESLAC (see Supplementary Information). 
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FIG. 4: Initialization and single-shot readout of a two nuclear spin qubit quantum register, (a), Experimental sequence. 
(b),(c), PL time trace recorded by continuously repeating the sequence with a magnetic field B — 5280 G applied along the 
NV defect axis. Each data point corresponds to approximatively 3000 repetitions of the readout sequence (20 ms). Quantum 
jumps linked to (b) the weakly coupled '^^C nuclear spin and to (c) the intrinsic ^''N nuclear spin of the defect can be 
distinguished, (d), Normalized nuclear-spin dependent photon counting distributions measured with the initialization thresholds 
set to A/i,darfc = 110 counts and A/i, bright ~ 160 counts. The projective readout fidelity is optimized for a discrimination threshold 
■N'r,opt = 135 counts. 



ing to depolarization at a rate 

1 ■ 

Considering this process at the main source of depolar- 
ization, the nuclear spin relaxation time might exhibit 
a quadratic dependence with the appUed magnetic field. 
The experimental results depicted in Figure [3] confirm 
this behavior at high fields, while two drops can be ob- 
served around B ~ 510 G and B ^ 1020 G, correspond- 
ing to level anti-crossings in the excited state and in the 
ground state, respectiveljJ22] (Fig. [ija)). Around such 
magnetic field strengths, the secular approximation is not 
valid and additional electron-nuclear spin flip-flop terms 
A±\S 1+ + S+I-]/2 need to be added to the Hamil- 
toniaiP^'^, where A± = {Axx + Ayy)/2. As shown 
in Fig. [3j the experimental data are well fitted by a 
simple model including depolarization induced by the 
anisotropic hyperfine interaction and spin mixing at the 
level anti-crossings (see Supplementary Information for 
details). We note that the bright state always exhibits a 
longer relaxation time than the dark state. Furthermore, 
this effect is independent on the nuclear spin state (|t) or 
ID) used as control state in the CNOT gate. When the 
dark state is detected, a shorter nuclear spin lifetime is 
always observed because in this case the system spends 
on average more time in the rug = — 1 electronic spin 
sublevel, for which the anisotropic component of the hy- 
perfine tensor induces nuclear spin flips. 

According to equation Q, a long nuclear spin life- 
time can be observed either for a ^^C nuclear spin with 
a weak anisotropic component of the hyperfine interac- 
tion, i.e. placed on a lattice site with a small angle with 
respect to the NV defect axis, or for an applied mag- 
netic field such that 7„_B ^ {Azz, Aani)- In a diamond 



sample with a natural abundance of ^'^C isotope (1.1%), 
the ESR linewidth is on the order of 200 kHz, which 
puts a limit to the weakest detectable hyperfine coupling 
strength in conventional ESR spectroscopy (Fig. 1(c)). 
Apart from the ^"^C nuclear spin studied in detail in this 
work, quantum jumps were also observed for a ^'^C cou- 
pling strength Azz — 380 ± 10 kHz (lattice site O in 
Ref. [28]) with a much shorter relaxation time (see Sup- 
plementary Information for details). For stronger hy- 
perfine coupling strengths, no quantum jumps could be 
observed for magnetic fields up to 5000 G. The proba- 
bility to find weakly coupled ^'^C nuclear spins would be 
significantly improved by using isotopically purified dia- 
mond samplei^^. However, we note that the speed of the 
single-shot readout measurement decreases with the ^^C 
coupling strength owing to the required spectral selectiv- 
ity of the quantum logic. 

Finally, we demonstrate single-shot readout in a two- 
qubit register by including the intrinsic ^"^N nuclear spin 
of the NV defect. For this experiment, the CNOT gate is 
performed with a single narrowband MW 7r-pulse which 
selectively drive the ESR transition for a given state of 
the register, e.g. state (Fig. 1(b) and Fig. 4(a)). 

The PL time trace then exhibits quantum jumps linked 
to both nuclear spin species, which can be easily dis- 
tinguished because their characteristic relaxation times 
differ by orders of magnitude (Figs. 4(b) and (c)). In- 
deed, although the ^*N nuclear spin shares its symmetry 
axis with the NV defect {Aani = 0), its relaxation time 
is only a few tens of milliseconds because a strong hyper- 
fine contact interaction in the NV defect excited-state 
A± « 40 MHz induces fast electron-nuclear spin flip- 
flop^. From the nuclear-spin dependent photon counting 
distributions, we infer that the two nuclear spin qubits 
can be initialized into state (dark state) with a fl- 
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delity higher than 98% by using an initialization thresh- 
old Ni^dark = 110 counts (Figs. 4(d)). We note that 
any state of the register could be deterministically pre- 
pared and readout by changing the frequency of the MW 
used for the CNOT gate. From the overlap between the 
photon counting distributions, we extract a projective 
readout fidelity = 83 ± 2%, limited by the ^"^N nuclear 
spin relaxation time. This value could be significantly 
improved by increasing the magnetic field strength in or- 
der to decouple more efficiently the ^^N nuclear spin from 
the electron spin dynamic^. 

The reported initialization and single-shot readout of 
two nuclear spin qubits combined with well-established 



techni ques o f coherent manipulation within the quantum 
registepl^'i^ pave the way towards tests of quantum 
correlations with solid-state single spins at room temper- 
ature^ and implementations of simple quantum error 
correction protocols^^. 
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SUPPLEMENTARY INFORMATION 



A. Experimental methods 

1. Experimental setup 

We study native NV defects hosted in a commercial [100]-oriented high-purity diamond crystal grown by chemical 
vapor deposition (Element6) with a natural abundance of ^^C isotopes (1.1%). Individual NV defects are optically 
isolated at room temperature using a confocal microscope. A laser operating at 532 nm wavelength is focused onto the 
diamond sample through a high numerical aperture oil-immersion microscope objective (Olympus, x60, NA=1.35) 
mounted on xyz-piezoelectric scanner (MCL, Nano-PDQ375). The red-shifted NV defect PL is collected by the 
same objective and spectrally filtered from the remaining excitation laser with a dichroic filter and a bandpass filter 
(Semrock, 697/75 BP). The collected PL is then directed through a 50-/im-diameter pinhole and focused onto a silicon 
avalanche photodiode (Perkin-Elmer, SPCM-AQR-14) operating in the single-photon counting regime. Laser pulses 
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are produced with an acousto-optical modulator (MT200-A0.5-VIS) with a characteristic rising time of 10 ns. For all 
experiments, the optical pumping power is set at 300 /iW, corresponding to the saturation power of the NV defect 
radiative transition. 



2. ESR spectroscopy 



The NV defect ground state has an electronic spin 5=1 that can be efhciently polarized into its iris = sublevel 
through optical pumping^. In addition, the PL intensity is significantly higher 30%) when the mj = state is 
populated allowing the detection of electron spin resonances (ESR) on a single NV defect by optical meanfl. 

Coherent manipulation of the NV defect electron spin is performed by applying a microwave field through a copper 
microwire directly spanned on the diamond surface. Electron spin resonance (ESR) spectroscopy is performed through 
repetitive excitation of the NV defect with a resonant microwave 7r-pulse followed by a 300-ns read-out laser pulse^. 
ESR spectra are recorded by continuously repeating this sequence while sweeping the 7r-pulse frequency and recording 
the PL intensity. The microwave power is adjusted in order to set the 7r-pulse duration to 4 /is, as verified by recording 
electron spin Rabi oscillations. In this conditions, the ESR linewidth is given by the inhomogeneous dephasing rate 
of the NV defect electron spin which is on the order of 200 kHz for a diamond sample with a natural abundance of 
"C isotope EI. 

The ESR spectrum of a single NV defect coupled with a nearby ^'^C nuclear spin shows six nuclear-spin conserving 
transitions (see Fig. 1(c) of the main text and Fig. [5](a)). Indeed, hyperfine interaction with the intrinsic -'^^N nuclear 
spin (/ = 1) leads to a splitting of = 2.16 MHz between ESR frequencies associated with different ^^N nuclear spin 
projections^. These lines are further split through hyperfine interaction with the "'^'^C nuclear spin (Fig. [HJa)). Even 
at high magnetic fields, the ^^N nuclear spin populations are unbalanced owing to dynamic nuclear spin polarization 
induced by optical pumpinj^ (see Fig. 1(c) of the main text and Fig. [sjja)). We note that this effect is responsible 
for the short polarization time of the ^'*N nuclear spin under optical illuminatiorP. 
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FIG. 5: Single-shot readout of a ^^C nuclear spin with a hyperfine coupling strength of 380 kHz. (a), Optically detected ESR 
spectrum recorded for a single NV defect coupled to a ^^C nuclear spin with a coupling strength Azz = 380 ± 10 kHz. The 
solid line is a fit with Gaussian functions, (b), PL time trace recorded by continuously repeating the readout sequence shown 
in Fig. 2(a) of the main text, with a magnetic field B = 2000 G applied along the NV defect axis. Each data point corresponds 
to approximatively 3000 repetitions of the readout sequence (20 ms). The solid line is a fit with a two states hidden Markov 
model, (c), Normalized nuclear-spin dependent photon counting distributions obtained with the initialization thresholds set to 
■N'i.i — 150 counts and A/i.t ~ 240 counts. With this values the initialization fidelity is 94%. 
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3. Magnetic field alignment 

A permanent magnet mounted on a xyz-translation stage is used to apply a static magnetic field along the NV defect 
axis. Preliminary alignment of the field is done by optimizing the PL intensity because any off-axis components of the 
magnetic field quench the NV defect PlI^Sl. The field alignment is then more precisely realized by measuring the sum 
E of the resonance frequencies and i^-i, linked to the transitions = ^ ms ~ +1 and = 0—7- mg = — 1, 
respectively. For a perfectly aligned magnetic field, S = v^i + v^i = 2D, where D is the zero-field splitting. In our 
experiments, this criteria is completed with a precision of about 100 kHz, corresponding to a magnetic field alignment 
with a precision better than 0.2° for a magnetic field oi B = 2000 G. 



B. Single-shot readout of a ^'^C nuclear spin with a hyperfine coupling strength of 380 kHz 

As indicated in the main text of the manuscript, quantum jumps were also observed for a ^^C coupling strength 
Azz — 380 ± 10 kHz, as shown in Figure [sjjb). For a magnetic field B = 2000 G applied along the NV axis, the 
characteristic relaxation times of the ^^C nuclear spin are Ti_^(^i,right) = 69±4 ms and Ti^^(^^ark) = 35±2 ms. From the 
overlap between the nuclear-spin dependent photon counting distributions (Fig. [sf^c)), we infer a projective readout 
fidelity J" = 77 ± 3%. 



C. Nuclear-spin dependent photon counting distributions 

In this section we describe how the histograms in Figure 2(c) of the main paper are obtained and how the statistics 
of the counted photons in our single-shot readout process is modeled. In the following the ^^C nuclear spin states 
are denoted as bright and dark, corresponding to a high and a low photon-counting signal, respectively. These states 
correspond to either |t) or ||) depending on the nuclear-spin conserving ESR transition used for the CNOT gate. 
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FIG. 6: (Left) Distribution of consecutive single-shot measurements. (Right) Histograms of photon counting distributions for 
dark and bnght states by using Ni^dark = 615 counts and Nbright ~ 845 counts. 



In order to obtain the nuclear-spin dependent photon counting histograms, we first consider the distribution of 
two consecutive measurements, Af{k) and Af{k + 1), from the PL time trace obtained by applying continuously the 
single-shot readout sequence. As it can be seen on Figure |6j this distribution is highly concentrated in two regions 
representing both nuclear spin states, dark and bright. We define the threshold value Mi. dark (resp. Afi^bright) on the 
measurement f\f{k) to declare that the nuclear spin state is initialized in the dark (resp. bright) state. Then, for 
all measurements that satisfy Af{k) < Mi, dark (resp. Af{k) > Mi, bright), a histogram is constructed from the set of 
points {M{k + 1)} (see right panel in Figure [6]). The lower (resp. larger) the Mi^dark (resp. Mi^bright) threshold, the 
smaller the initialization error but the larger the uncertainty as we obtain fewer points to built the histograms. As a 
good compromise between these two effects, we choose Mi, dark = 615 counts and Mbright = 845 counts. Finally, each 
distribution is binned in intervals of 15 counts to construct the nuclear-spin dependent photon counting histograms 
shown in Fig. 2(c) in the main paper. 
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We now explain how these histograms are modeled. As described in section D of this supplementary information, 
the flipping rate of the nuclear spin depends mainly on the anisotropic part of the hyperfine interaction between the 
electronic spin and the ^^C nuclear spin. This causes an increase of the flip-flop rate when the electronic spin of the 
NV center is in state rris — ±1. As a consequence, the dark state always presents a smaller relaxation time Ti than 
the bright state when the readout sequence is applied continuously. 

Therefore we model the flip-flop events with a two-rate poissonian distribution. The time Yd (resp. Yg) the nuclear 
spin spends on the dark (resp. bright) state distributes exponential with rate Xu (resp. A^). Following a similar 
procedure to that described on ref.H^, we model the number of photons detected over a measurement time T with a 
random variable Z, by considering the statistics of the photons associated with each nuclear spin state when there is 
no flip, one flip and two flips over the measurement time. 

When there is no-flip, the number of photons is given by the random variable, 

Zi = Xi , 1 = {B, D} (5) 

where Xj is a random variable that distributes normal, Ni{Xi — Xi), with mean /i/ — (resp. fii = f^s) and 
variance aj = cr|j (resp. aj = a^) if the nuclear spin state is dark (resp. bright). The distribution of the number of 
photons when there is no-flip /^^ (z) can then be calculated by taking the derivative with respect to z of the cumulative 
distribution, 

fl {z) = ^ I Ni{x,)P{Yj > T)dx, = Ni{z)e-^'^. (6) 

Jxi<z 

When there is one flip during the measurement time T, we model the number of photons by the random variable, 

Zj = ^Xj + ^r^Xj,, (7) 

where Yj distributes exponentially with constant A/. The set of indexes {/, /'} denotes the state of the nuclear spin 
and can be either {D, B} or {B, D}. The distribution is given by 

fzjiz) = ^ J Ni{xi)Np{xi,)P{Yi = t)P{Yi, >T- t)dxjdxi,dt (8) 

^xi + [l— ^)x <z 



H.. [\ A/Te-^^^V-^^'^d-") r [z-{u^lI + {l-u)^JiI 



^2^[u^aj H- (1 - u)2a|,] I 2[u^<yj + (1 - u)^^!, 



(9) 



Note that when A/ = A// = A, the distribution becomes /^(z) = Pj duN{fi, a^), where Pf = XTe iV(^, a^) is 
the Normal distribution with mean n and variance cr^, /i — u^j + (f — u)^j' and = a'jv? + (1 — uf'a'j,. 
Similarly, we model the two-flip case by a random variable 

.|x, + l£^x„ + I^x, ,10) 
i_!£_2i),, + ii^A> ,11) 

with distribution 

,2 , ^ [\ \iT\i,T{l-u)e-^'^e+<^'^->^''^) \ [z - {u^iI + {I - ^l)^lI,)] 



^2ti[u^(j] + (1 - m)V2,] I 2[uV2 + (f-7i)2a2, 



(12) 



Note that when A/ = A// = A, the distribution /|(z) = Pf du2{l - u)N{fi,a'^), where Pf ^ (AT)2/2e-^^, 
/X = ufii -I- (1 — u)fiii and cr^ — ,72^2 _(- (1 _ u)^aj,. 

Finally, we fit the experimental photon counting distributions shown in Fig 2(c) of the main paper to 

Siiz)^fl{z) + fhAz) + .fl{z) , (13) 
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where the fitting parameters are the average photon counting numbers ijd and /zb, the flip-flop probabiHties TXd 
and TXb for the nuclear spin in state dark and bright, respectively. On t he other hand, it is known that three level 
systems with a metastable state exhibit super-poissonian charactei^^ESl. As the photons associated with the dark 
distribution involve the passage of the electron throught the metastable singlet state ^Ai, the variance of the dark 
distribution is left as a fitting parameter, meanwhile the variance of the bright distribution is set to = ^b- As 
a result we obtain ~ 637 counts, — 816 counts, TA_d = 0.084 and TXb = 0.0485, leading to Ti^dark — 1-43 s 
and Ti^iij-igfif = 2.47 s in fair agreement with our experimental results (Ti^^j.^, — 1.5 i 0.2 s and 7i f,j,jg/j( = 2.4 ±0.1 
s). ' 



D. Nuclear-spin depolarization processes 

In this section, we discuss the depolarization processes of the ^'^C nuclear spin during the multiple repetition of the 
single-shot readout sequence. 



1. System Hamiltonian 

The system consists of a single NV defect coupled by hyperfine interaction with its intrinsic ^"^N nuclear spin and a 
nearby "'^^C nuclear spin. Since we focus on the evolution of the ^"^C nucleus, we do not consider the interaction terms 
linked to the ^^N. Assuming a magnetic field B perfectly aligned along the NV defect axis, denoted as the z-axis, the 
system Hamiltonian, in both the ground and excited states, reads 

=H«+H«-fHil (14) 
with the index i refcring cither to the ground state {gs) or to the excited state (es), and 

< n'i^ =inBi, (15) 
^ uf_^ = s • ^(') • i , 

where is the zero-field splitting of the NV defect electronic spin - D^'^") ~ 1.42 GH^^^and D'^as) ^ 2.87 GHz - , 
7e ~ 2.80 MHz.G^^ and 7„ ~ 1.07 kHz.G"^ are respectively the gyromagnetic ratio of the electronic spin and of the 
^^C nuclear spin, and A^"^^ is the hyperfine tensor. 

Given the small value of the hyperfine coupling strength considered in the main text of the manuscript {A\\ ~ 258 
kHz), the ^^C must be located few lattice sites away from the NV defeclP^. However, no correspondence with a 
specific lattice site of the diamond matrix is available owing to the current accuracy of ah initio calculations of the 
electronic spin wave functiorP. In the following, we neglect the contact term of the hyperfine interaction and thus 
assume a purely point-like dipolar interaction. The hyperfine tensor is therefore considered identical in the ground and 
in the excited states y^'-'^ — A. Furthermore, the significant ^^C polarization at the excited state level anti-crossing 
suggests the polar angle between the ^^C lattice site and the NV defect axis is smalP^. This is further supported by 
the observed long polarization time of the ^^C nuclear spin while applying repetitive readout. Thus, we can neglect 
the terms of the hyperfine interaction proportional to S'±/±'^and the Hamiltonian can be approximated as 

n^'^ ~ Hi')+Hi')+A.^.4 + ^[^+/.e-'^ + 5_4e+*^] (16) 

+ r^[S+i_ + SJ+] + ^[Sj+e-''^ + Sj-e+''-^] (17) 

where Aani = \jA?^.j. + A?^y , tan0 = Azy/Azx and Aj_ = {Axx + Ayy)/2 = —Azzl'^ because the hyperfine tensor is 

traceless for a pure dipolar interaction. The last two terms of this Hamiltonian are responsible for the undesired ^^C 
nuclear spin flips. 



2. Nuclear spin flip processes 



The first nuclear spin-flip term is proportional to the perpendicular component of the hyperfine tensor A±, and 
connects the state | — l)e It) [resp. | 0)e |t)] with the state | 0)e \l) [resp. | -fl)e as depicted on Figure[7]ja). The 
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FIG. 7: Schemes illustrating the coupling between the different states | ms)e |ti)i resulting either from the perpendicular 
component A± (a) or from the anisotropic component (b) of the hyperfine interaction, (c), Nuclear spin-flip rates linked to 
the anisotropic hyperfine coupling (red curves), and to the perpendicular hyperfine coupling in the ground (green curve) and 
excited state (blue curve). The parameters are set to Azz = 258 kHz, 7„ = 1.07 kHz.G"\ je = 2.80 MHz.G"\ D'*"' = 2.87 
GHz and D^""'' = 1.42 GHz and Aani = 1 kHz (red solid line), 10 kHz (red dashed hne) and 100 kHz (red dotted fine). 



nuclear spin-flip rate induced by this coupling is proportional to the transition probability between the two interacting 
states. Since the detuning between | 0)e |t) and | +l)e \i) is always larger than the one between | — l)e It) and | 0)e |4,), 
the main depolarization rate linked to the perpendicular component of the hyperfine tensor can be expressed as 

2.4^ 

(18) 



1± = 



1 

r>c 



2A\ + (D^^ - jeB) 



As shown on Figure [tJ^c) , this rate becomes significant when the system gets close to level anti-crossings, i.e. for 
jgB, occurring around B ^ 1020 G and B ^ 510 G, in the ground state and in the excited state, respectively. 



The second nuclear spin- flip term in equation (17) is proportional to the anisotropic component of the hyperfine 
tensor Aani (Fig. [rjjb)). For all the experiments reported in the main text, we use nuclear-spin conserving transitions 
between the rus = and = — 1 electron spin manifolds. Since the anisotropic hyperfine interaction couples the 
states I — l)e It) and | — l)e the nuclear-spin-flip rate is therefore given by 



1 



^2 . 

am 



Ti., 



(19) 



A^ 

(A. - 7nS)2 

The evolution of this rate with the magnetic field amplitude B is shown on Figure [7]jc) for different values of Aani- 
Even for a value as small as Aani = 1 kHz, the anisotropic component of the hyperfine interaction is always the 
dominant depolarization process, except near the two level anti-crossings. This is supported by the experimental 
results depicted in Figure 3 of the main text, which indicate that the bright state, corresponding to the electronic spin 
in state | 0)e, always exhibits a longer relaxation time than the dark state. 



3. Ti evolution with the magnetic field 



The evolution of the nuclear spin relaxation time versus the magnetic field is depicted in Figure 3 of the main text. 
At a given magnetic field amplitude, Ti^^(^t,right) (red dots) and Ti^^i^^ark) (blue squared dots) are inferred from a fit 
to a PL time trace showing quantum jumps with a two-state Hidden Markov ModeP*^. 

We simply consider the total depolarization rate 71 as a weighted average of the three depolarization rates ^ani, 
leading to the formula 

1 



and 7_l 



71 



a, 



am Jam \ /_|_ 



(es) (es) 



(20) 
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where {aani,cif^\a^^'^^} are coefBcients linked to the optical pumping power, the intrinsic photophysical parameters 
of the NV defect and the time the system spends in each state | i)e \ti) during the single-shot readout sequence. The 
solid lines in Figure 3 represent a fit of the experimental results using equation ( 20 1 with {aani , C(±^^ , Q^Il } as fitting 
parameters and setting Aani = 10 kHz. Although this highly simplified model does not allow to extract quantitative 
information, it reproduces fairly the general trend of the experimental data. The development of a more precise model 
would require to introduce the NV defect dynamics under optical pumping, including ionization of the defect in the 
neutral charge state NV" (Ref.Sil). 

We note that the spin relaxation time could be significantly enhanced by aligning the magnetic field along the ^^C 
hyperfine field rather than the NV axis, leading to Aani ~ 0. This could be realized for a ^'^C placed at a lattice site 
with a small polar angle with respect to the NV defect axis in order to avoid any significant electronic spin mixing 
which degrade the ESR contrast. 
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